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We have developed a theory of exciton-polariton band structure of resonant three-dimensional 
photonic crystals for arbitrary dielectric contrast and effective mass of the exciton that is excited 
in one of the compositional materials. The calculation has been carried out for a periodic array of 
spheres embedded in a dielectric matrix. It has been shown that the position of the lower branches of 
the polariton dispersion curve monotonously depends on the exciton effective mass and is determined 
by the coupling of light with the first few states of the mechanical exciton quantum confined inside 
each sphere. Particularly, we have studied the role of excitonic effects on the photonic-crystal band 
gap along the [001] direction of the Brillouin zone and presented an analytic description of the 
polariton dispersion in terms of the two-wave approximation. 



PACS numbers: 71.35.-y, 71.36.+C, 42.70.Qs 
1. INTRODUCTION 



The concept of photonic crystals was put forward by 
Yablonovichi and Johri^ in 1987. Since then this term is 
used for media with the dielectric susceptibility varying 
periodically in the space and allowing the Bragg diffrac- 
tion of the light. The theory of photonic crystals has 
been deve loped in a number of subsequent works, see, 
e.g., Refs. 3 4 5 6 7. The main goal of these studies was 
to determine the photon band structure and analyze the 
sequence of allowed bands and stop-bands (gaps) for dif- 
ferent directions of the wave vector in the first Brillouin 
zone. The simplest model realization of a photonic crys- 
tal is a structure grown from two materials, A and B, with 
different dielectric constants ea and Eb'- a periodic lay- 
ered medium ...A/B/A/B... in case of one- dimensional 
photonic crystals and periodic arrays of cylinders and 
spheres of the material A embedded in a dielectric ma- 
trix B, in case of two-dimensional and three-dimensional 
photonic crystals, respectively. Periodic structures where 
the dielectric susceptibility of one of the compositional 
materials, as a function of the frequency lo, has a pole at 
a certain resonance frequency are grouped into a special 
class of resonant photonic crystals. In such systems the 
normal light waves are excitonic polaritons. In Refs. 00 
the dispersion of light waves has been calculated tak- 
ing into consideration the frequency dependence of the 
dielectric susceptibility in the frame of the local mate- 
rial relation D = ea{<->j)E between the electric displace- 
ment and the electric field. In Ref. ^] the dispersion of 
excitonic polaritons in a resonant photonic crystal has 
been calculated taking into account only one level of the 
quantum-confined exciton in a sphere A and neglecting 
the difference between the dielectric constant e b and the 
background dielectric constant e a of the material A. In 
the present work we have studied theoretically the disper- 
sion of excitonic polaritons in a resonant photonic crys- 
tal making allowance for all available exciton quantum- 
confinement levels as well as for the dielectric contrast, 
i.e., for e a ^ e B - 



2. PROBLEM DEFINITION AND METHOD OF 
CALCULATION 

In this work the theory of resonant photonic crystals 
is developed for a periodic array of spheres A arranged 
in a face-centered cubic (FCC) lattice and embedded in 
the matrix B. The structure under consideration is char- 
acterized by seven parameters: R, a, Eb, E a , wur, w o 
and M. Here R is the radius of the spheres A, a is the 
lattice constant for the FCC lattice, Eb is the dielec- 
tric constant of the matrix, loq, wlt and M are the res- 
onance frequency, longitudinal-transverse splitting and 
translational effective mass of the triplet Is exciton ex- 
cited inside the spheres A, e a is the background dielectric 
constant that contains contributions to the dielectric re- 
sponse from all other electron-hole excitations. Thus, we 
assume the dielectric function of the bulk material A to 
have the form 

e a (uj, q) = £ a H — , u) CKC {q) = uj q + — (1) 

u} cxc (q)-LU 2M 

and, in that way, take into account both the frequency 
and spatial dispersion, i.e. the dependence on the light 
frequency u and wave vector q. The radius R is chosen 
so that, on the one hand, the spheres A do not overlap, 
i.e. R < a/2V2, but, on the other hand, R should exceed 
the Bohr radius of the ls-exciton in the material A and, 
hence, the exciton can be considered as a single particle 
with the mass M. In the following we ignore the fre- 
quency dependence of the parameters Eb and e a . More- 
over, hereafter we assume the material A to be isotropic 
and take into consideration the bulk Is exciton states 
only. It follows then that the problem is reduced to the 
solution of a system of two vector equations, namely, the 
wave equation 

rotrot£;(r) = [e(r)E{r) + 4nP cxc (r)] (2) 

and the material equation for the ls-exciton contribution 
to the dielectric polarization 

(-i A + ^-)^cW = ^W, (3) 
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where e(r) — e a inside the spheres and e(r) = eb out- 
side the spheres, E(r) and P exc (r) are the electric field 
and the excitonic polarization at the frequency lo. On a 
spherical surface separating materials A and B we impose 
the standard Maxwell boundary conditions: continuity of 
the tangential components of the electric and magnetic 
fields, and the Pekar additional boundary condition for 
the excitonic polarization: the vanishing vector P QXC (f) 
at |r — a | = R, where the transnational vectors a define 
the centers of spheres A. 

Due to the periodicity of the structure we can seek 
solutions of Eqs. (J2J and © in the Bloch form satisfying 
the condition 



E k (r + a) = c ika E k (r). 



(4) 



Here k is the exciton-polariton wave vector defined 
within the first Brillouin zone. We remind that, for a 
FCC lattice, the latter is a dodecahedron bounded by six 
squares and eight hexagons. 

Below we present the results of calculation of the 
exciton-polariton dispersion aj nk , where n is the branch 
index. The computation was mainly performed by us- 
ing a photon analogue of the Korringa-Kohn-Rostoker 
(KKR) method 11 ! 12 ^ 3 . In this method (i) the electric 
field is decomposed in the spherical waves, or more pre- 
cisely, in the vector spherical functions centered at the 
points r — a, and (ii) following the consideration of the 
light scattering by a single sphere and the introduction of 
a structural factor the dispersion equation is transformed 
to 



G j ' m > (T ' < j m(T (k,u;)Rj (7 (Lj)\ = . (5) 



Here Rj a are the coefficients describing the scattering of 
spherical waves by a single sphere A, they depend on the 
total angular momentum j and the polarization index a 
discriminating the magnetic and electric spherical har- 
monics but are independent of the angular-momentum 
component m. Note that, for a spherical scatterer, these 
coefficients relate the incident field Eo(r) cx Jj mo -(r) 
with the scattered field E sc (r) oc Hj m<J (r), where 
Jjma-, Hj ma are the vector spherical functions^ 3 -. Ajiki 
et alm^ have calculated values of Rj a (u) taking into 
account the exciton resonance in the spheres A and 
the finite exciton effective mass M. In contrast to the 
scattering matrix for a single sphere which is diagonal, 



Rj a 5 



•j'm'a' .jma, the matrix of structural 
factors Gj'm'a-' j ma (k,uj) has nonvanishing diagonal and 
off-diagonal components. It should be mentioned that 
the both matrices are frequency dependent whereas only 
the structural matrix G depends upon the polariton wave 
vector k. At the same time G is independent of the ex- 
citonic para meters and coincides with the matrix consid- 
ered in Refs. lllliail where the exciton states are disre- 
garded. 

In addition to the KKR method, in section 4 we apply 
the Green-function technique for the analysis of separate 
contributions of quantum-confined excitonic levels to the 




FIG. 1: Exciton-polariton band structure of a photonic crys- 
tal with a FCC lattice of spheres A inserted into the matrix B. 
The calculation is performed neglecting the dielectric contrast 
and for the set of parameters indicated in the text. 



polariton dispersion, and in section 5 we use the two- wave 
approximation which allows an analytic description. 



3. POLARITON DISPERSION FOR A FINITE 
EXCITON MASS 

First we focus the attention upon purely excitonic ef- 
fects and neglect the dielectric contrast assuming e a — 
eb- Then in the absence of exciton-photon interaction, 
i.e. for ojlt = 0, the medium becomes optically homoge- 
neous and photons propagating therein are characterized 
by the linear dispersion to — cq/ns with the refractive in- 
dex riB = \f£B- In the reduced zone scheme, this single- 
valued relation between the frequency and wave vector q 
turns into the many-valued (or multi-branch) dispersion 



\k + b\/n B 



(6) 



where b is the reciprocal lattice vector reducing q to the 
vector k = q — b that lies in the first Brillouin zone. A 
nonvanishing value of wlt leads to a mixing between pho- 
tonic and excitonic states and formation of hybrid polari- 
ton excitations characterized by a complicated multiband 
dispersion uj nk . As a result the wave {n,k) is a mixture 
of two or more photonic states © with the same k but 
different b. 

Fig. 1 presents the dispersion of exciton-polaritons cal- 
culated for a FCC lattice and the following set of struc- 
ture parameters: 

£ a = £b = 10, R = a/4, fiUJi — 2eV , cj lt = 5 x 10~ 4 tui , 



•y/371 



1.1 , M = 0.5 m . 



\ LOinsa I 

Here mo is the free electron mass in vacuum and, instead 



3 



of the bulk exciton resonance frequency loq, we intro- 
duced the resonance frequency 



M — > but loi = const, this branch is remarkably af- 
fected by the lower level J2J only. For M — > oo, the 



2M 



(7) 



of the ground state of the exciton quantum-confined in 
a sphere of the radius R. Note that, for the structure 
tuned to the Bragg resonance u>% = cUl /rib at the point 
L of the Brillouin zone with k^ = ^/3n/a or u>\ — ckx /rib 
at the point X with kx — 2ir/a, the parameter P equals 
to 1 and 3v3/8 ~ 0.65, respectively. For P = 1.1, the 
anticrossing between the horizontal line u> = lo\ ( "bare" 
exciton branch) and the line u> = ck/ng ("bare" photon 
branch) occurs inside the Brillouin zone at k amounting 
approximately 97% of k^ and 84% of kx- For the sake 
of completeness the dependence u>(k) is shown not only 
along the high-symmetry crystallographic directions k \\ 
[001] (points A) and k || [111] (points A) but also along 
the straight lines X - W, W - K, X - U and U - L. 

In Fig. 1 the frequency region is cut from above in or- 
der to keep only few branches of the dispersion curve. 
In the cut region the dispersion is presented by a dense 
network of polariton branches which appear as a result 
of the anticrossing of "bare" photonic branches © with 
a set of close- lying exciton quantum-confinement levels. 
The network has a tangled character complicated for de- 
piction. Therefore, we concentrate here on the analysis 
of the role played by the photonic-crystal parameters in 
the formation of the lower polariton branches. 

In Fig. 2 the solid curves show the same disper- 
sion branches as in the previous figure but in an in- 
creased scale and in the vicinity of the points X and L. 
For comparison we also present, by dashed-and-dotted 
lines, the lower branch of the dispersion curve calculated 
taking into account only the ground exciton quantum- 
confinement level. The calculation was performed by 
reducing the value of M to 0.01 mo and decreasing ojq 
so that the frequency oji in Eq. Q was kept invariant. 
The dashed-and-dotted lines coincide with those calcu- 
lated by the method proposed for exciton-polaritons in 
Ref. ha Dashed curves illustrate the opposite limiting 
case of quite heavy excitons, M — > oo. In this case the re- 
lation between the exciton polarization and electric field 
becomes local: 



XE, 



X 



4-7T LOq 



It means that values of k corresponding to a certain fre- 
quency u can be found in the same way as it is done for 
nonresonant photonic crystal with the dielectric suscep- 
tibilities Eb and ea = E a + ^X- The calculation shows 
that the dashed line in Fig. 2 is practically indistinguish- 
able from the lower branch obtained according to Eq. © 
for M = 5mo. The lower polariton branch is formed as a 
result of "repulsion" of the photon branch with 6 = 
towards the long-wavelength side because of the interac- 
tion with the exciton quantum-confinement levels. For 





FIG. 2: Dispersion of exciton polaritons in a resonant pho- 
tonic crystal in the spectral region adjoining the resonance 
frequency of the lower exciton level u\ and for the wave vec- 
tors k || [111] (points A) and k || [001] (points A). The 
dispersion curves are calculated for the exciton effective mass 
M — 0.5 mo (solid lines), M — * oo (dashed lines) and M — > 
(dashed-and-dotted lines). Other parameters are the same as 
in Fig. 1. 



other levels act upon this branch to the maximum since 
in this limit their resonance frequencies coincide and are 
equal to luq. We conclude then that the lower polariton 
branch corresponding to the finite mass M must always 
lie between the dashed-and-dotted and dashed lines, in 
agreement with the results of calculation shown in Fig. 2. 



4. COMPARATIVE CONTRIBUTION OF 
INDIVIDUAL QUANTUM-CONFINED EXCITON 
LEVELS 



In the following we use the notation ujx (M) for the fre- 
quency at the X point in the lowest polariton branch for 
the exciton effective mass M in the material A, u)x(0) 
and ojx(oo) for values of this frequency for M — > 
and M — > oo, and, finally, uix for the arithmetic mean 
[wx(0) + uJx(oo)]/2. One can see from Fig. 2 that 
the difference oj\ — lux noticeably exceeds the difference 
k>x(0) — wjt(oo). This means that the main contribu- 
tion to the position of the frequency lux (M) should come 
from the exciton-photon coupling with the ground level 
0. In this section we analyze in detail the influence of 
the ground and excited levels of the quantum-confined 
exciton on the lower polariton branch. To this end we 
follow Ref. El use the Green-function technique and 
expand a solution of equation with arbitrary function 
E(r) in a series 
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P«c(r) = |i£;-^ -a) / K(r'-a)E(r')d 3 r' (8) 

r'-a|<i? 



over the eigen states of the mechanical exciton in a 
sphere with an infinitely high barrier. Here the index 
v = (n r ,l,m) characterizes the exciton state and com- 
prises, respectively, the radial quantum number, angular 
momentum and its projection on the z axis, w„ is the 
exciton resonance frequency in the state v. It is worth- 
while to note that the functions Q v (r) satisfy Eq. with 
uj = lo v and without the inhomogeneous term, i.e., with 
E{r) = 0, they are normalized to unity and correspond to 
the sphere centered at the point r = 0. The frequencies 
w„ are conveniently presented as ujo + (h/2MR 2 )x^ lf 
where i„ ri i are dimensionless numbers. For the few low- 
est energy levels, their values are given by tt (Is), 4.493 
(lp), 5.764 {Id), 2tt (2s), 6.988 (1/), see e.g. Ref. E 
with symbols in parenthesis indicating the radial num- 
ber n r and the orbital momentum, s for I — 0, p for 
I = 1 etc. Note that these symbols characterize the exci- 
ton quantum-confined state whereas the internal exciton 
state is Is all the time. Substituting the expansion © 
into the wave equation J2J and keeping in the sum over 
v any one term or a finite number of terms we can study 
the effect of these terms on the formation of the polari- 
ton dispersion. The results are illustrated in Fig. 3. It 
is seen that (a) the position of the frequency u>x(M) is 
mainly determined by the exciton level Is, (b) the po- 
lariton frequency at the X and L points is unaffected 
by the lp level, and (c) it suffices to add the contribu- 
tions from the lei and 2s levels in order to describe quite 
well the position of uix(M) whereas, at the L point, the 
sum converges slowly. Obviously, this difference in the 
convergence is related to the choice of P = 1.1 in which 
case the effect of anticrossing near the L point is much 
stronger compared with that at the X point. 



5. POLARITONS IN PHOTONIC CRYSTALS 
WITH A DIELECTRIC CONTRAST 

The dashed lines 1 and 2 in Fig. 4 depict the disper- 
sion of light waves in the vicinity of the X point of the 
Brillouin zone in a nonresonant photon crystal. Exactly 
at the X point these waves are characterized by the sym- 
metry X$ and X$t. The solid lines 1', 2' and 3' show 
the dispersion branches calculated taking into account 
only the exciton-photon coupling with one excitonic level 
Is and choosing its resonance frequency in the middle 
between the frequencies of the "bare" photons X$ and 
X51 . An additional branch 4' represents longitudinal ex- 
citon states which are practically dispersionless, in the 
following it is not discussed. Finally, the dashed-and- 
dotted lines 1", 2" illustrate the exciton-polariton dis- 






FIG. 3: Dependence of the exciton-polariton spectrum on the 
number of excitonic levels taken into account while calculating 
the dispersion curve by the Green-function technique. Lines 
1, 2, 3 and 4 are obtained with allowance made, respectively, 
for one, two, three and four lower levels; line 5 represents the 
exact result taking into account all levels of exciton quantum 
confinement. The structure parameters used are the same as 
in Fig. 1. 



persion with allowance made for all exciton levels. For 
the illustration we chose a photonic crystal with com- 
paratively weak dielectric contrast, eb = 12, e a = 13, 
in order to have the lowest stop-band in the direction 
k || [001], or the splitting between the X$ and X& states, 
comparable with the matrix element of exciton-photon 
interaction. It follows from Fig. 4 that the allowance for 
the Is quantum-confined exciton states optically active 
in the polarization E _L z results in a replacement of two 
branches 1, 2 by three branches 1', 2' and 3'. The ad- 
dition of contributions due to other (excited) excitonic 
levels leads to a transformation of the branch 1' into 1" 
and an appearance an extra branch 2". Simultaneously 
a dense network of polariton branches is formed in the 
frequency region uj > loq. It is not shown in the figure 
in order to avoid overloading the image. Instead, only 
the exciton-polariton eigen frequencies w v > ujq at the X 
point are indicated by short horizontal lines intersecting 
the vertical line X. 

In what follows we propose a two-wave approximate 
description of the polariton spectrum which is valid for 
a weak dielectric contrast and allows to understand the 
nature of branches indicated in Fig. 4 by primed and 
double-primed integers. Earlier this approximation was 
used in the analysis of the photon dispersion in nonreso- 
nant photonic crystals^. 

The Bloch solution 10} can be expanded over the space 
harmonics with the wave vectors k + b. For the vec- 
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[(fc 2 /fc ) 2 - e] £ 2 = £ '^i + E ^ E 4 a) *^n 
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The notation used are as follows: ko = u)/c, 

i = — [ e(r)d 3 r , e' = — [ e 4lTiz/a e(r)d 3 r , 
Vq J v J 



\/vo J 



£«WLT 



(11) 



FIG. 4: Photon band structure for a crystal with the dielectric 
contrast, e a = 12, eg = 13. Dashed lines are calculated for a 
photonic crystal neglecting excitons; solid lines are obtained 
taking into account only the lowest excitonic level; dashed- 
and-dotted lines and short straight lines are obtained taking 
into account all quantum-confined states of the mechanical 
exciton. 



tors k || [001] lying close to the X point, we keep 
in this expansion only two terms with k\ = (0,0, fci), 
k 2 = (0, 0, k 2 ) and k\ — k 2 — 4ir /a. At the X point 
we have k\ = — k 2 = 2ir/a. Thus, the electric field is 
approximated by 



E(r) = £; ie ifclZ + E 2 e 



(9) 



For doubly-degenerate polariton states of the A5 sym- 
metry compatible with the representations X$ and X51, 
the amplitudes E\ and E 2 are parallel to each other and 
perpendicular to the axis z | [001]. Substituting Eqs. © 
and into the wave equation , multiplying its terms 
by exp(— ikjz) and integrating over the unit primitive 
cell vq = a 3 /4 of the FCC lattice we obtain 

[(h/ko) 2 - e] Eh. = e'E 2 + E ^ E 1™) 



and we assume the origin of the Cartesian coordinate 
system to be chosen in the center of one of the spheres 
A. 

In the absence of the exciton-photon coupling and di- 
electric contrast, the right-hand sides of equations i|10|) 
vanish, e = Eb, and we arrive at two branches of the 
dispersion relation © for "bare" photons. At the X 
point these branches converge. In the presence of dielec- 
tric contrast, e' 7^ and the fourfold degenerate state 
at the X point splits into doubly-degenerate states X& 
(Ei = E 2 ) and X 5 (E\ = —E 2 ) with their eigen frequen- 
cies being 



, v , ck x , , ck x 
ye - s y. 



£ + £' 



(12) 



and the splitting being w(X$) —uj(Xsi) w (e'/e)(ckx /n), 
where kx — In /a, n — ye. The approximate equations 
(|12fl reproduce with high accuracy results of the exact 
calculation presented in Fig. 4 by dashed lines. 

In the absence of dielectric contrast, with allowance 
for the exciton level Is only and negligible difference be- 
tween the integrals and l[ 2 J the set of equations 110|) 
reduces to 



[(fci/fco) 2 - eg] E x = [{k 2 /k ) 2 - eg] E 2 = T x Jl s {E x + E 2 ) 



(13) 



r 



where I\ s is the integral (|ll(l corresponding to the wave 
vector kx- Provided that the frequencies u)\ and u>x = 
ckx/ng are close to each other, the photonic states A5, 
X51 are replaced by three doubly-degenerate polaritonic 
states: one of the symmetry X 5 with the frequency lux 
and two of the symmetry X51 with the frequencies 



u> = ± 



UJi - UJ X 



5 2 



(14) 



where S = y/uiitJurlis . For coinciding frequencies u>\ 
and u>x, the polariton spectrum in the direction k | [001] 



has a stop-band of the width 25 centered at u> = co±. The 
detuning of u>x from u>\ leads to shifts of the stop-band 
edges according to Eq. and, in addition, an allowed 
band is formed in the center of the stop-band, similarly as 
it happens in the resonant Bragg quantum-well structure 
(seeii and references therein) . 

For an approximate description of the solid lines in 
Fig. 4 it suffices to keep in the sums over v in Eqs. 1|10|) 
only the contribution due to the exciton quantum- 
confined ground state. In this approximation the mixing 
of the photon and exciton states of the symmetry X$i 
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leads to their repulsion and formation of hybrid waves 
with the frequencies 



± 



LUl - Lu(X 5 i) 



5 2 



(15) 



The previous equation fTljl is a particular case of this 
more general equation. Since the Is quantum-confined 
exciton does not interact with the light wave of the sym- 
metry X 5 , the X 5 photon frequency remains unchanged 
in the considered approximation. This allows to under- 
stand the closeness of the J^-polariton frequencies cal- 
culated neglecting the exciton effects and taking into ac- 
count only the Is level, see the points of intersection of 
the lines 2 and 2' with the vertical line X. Allowance 
for the excited excitonic levels in Eq. I|1U|) results in a 
shift of the frequency of the lowest polariton branch AV 
downwards (line 1"). However, their influence is small 
in comparison with the Is exciton. At the same time, in 
the formation of the polariton X5 (line 2") the main role 
is played by coupling of the photon A5 with the exciton 
lp and m — 0. The calculation shows that the branch 2" 
is well described by the two-wave model IjlOl) taking into 
account only one excitonic state (lp, m = 0). 



6. CONCLUSION 

We have developed a theory of resonant three- 
dimensional photonic crystals made up of two compo- 
sitional materials A (spheres) and B (dielectric matrix) 
for an arbitrary value of the exciton effective mass M and 
arbitrary dielectric contrast determined by the difference 
between the dielectric susceptibility Eb of the matrix and 
the background dielectric constant e a of the material A. 



For a finite value of M, the lowest polariton branch 
lies between the branches calculated in the two limiting 
cases, namely, in the absence of space dispersion, i.e., 
for M — > 00, and with allowance for only one excitonic 
level, i.e., for M — > and uj\ = const. For a satisfactory 
description of the polariton branches in the frequency 
region u> < wi, it suffices to take into account the light 
coupling with few lowest states of the mechanical exciton. 

For a resonant photonic crystal with the dielectric con- 
trast, we have analyzed the exciton-induced modification 
of the stop-band in the [001] direction if the exciton res- 
onance frequency u)\ is chosen to lie in the middle of 
the stop-band of the analogous photonic crystal with- 
out excitons. If the dielectric contrast is weak so that 
\sb — £a\ *C Eb then one can use the two- wave approx- 
imation that allows describe with good accuracy the re- 
sults of numerical calculation. The main contribution 
into the position of the lower edge of the polariton stop- 
band along the T — X line comes from the exciton Is 
states interacting with the light wave of the X51 symme- 
try. The upper edge of the stop-band is governed by the 
coupling of the exciton (lp,m = 0) with the photon of 
the X 5 symmetry. 

It should be noted that the developed theory with 
M — > 00 can be applied for the study of photonic 
crystals infiltrated with dye aggregates periodically dis- 
tributed in the space and characterized by a certain res- 
onance frequency of optical transitions, see— and refer- 
ences therein. The theory can be applied as well for the 
calculation of the dispersion of exciton-polaritons in res- 
onant two-dimensional photonic crystals, e.g., in periodic 
arrays of cylinders A embedded into the matrix B. 

The work is supported by the Russian Foundation for 
Basic Research (grant 05-02-16372) and the program of 
the Russian Ministry of Science and Education. 
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